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Abstract 

Calculations are carried out for the energy accommodation coefficient at a gas-surface interface 
using a recently developed classical mechanical theory of atom-surface collisions that includes both 
direct scattering and trapping-desorption processes in the physisorption well of the interaction 
potential. Full three-dimensional calculations are compared with the available data for the accom- 
modation of rare gases at a tungsten surface and good agreement is found for the heavier gases for 
which classical physics is expected to be valid at all measured temperatures. 
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I. INTRODUCTION 



The exchange of energy between a gas in contact with a surface is often characterized in 
terms of the energy accommodation coefficient. Although the origins of the concept of an 
energy accommodation coefficient can be traced to J. C. MaxwelU it is Knudsen who gave 
it a proper physical definition under what are now known as the conditions of rarefied gas 
dynamics.— 

The Knudsen energy accommodation coefficient has values that range from zero to unity, 
with a value of unity arising if the gas achieves equilibrium with the surface after colliding 
with it and a value of zero implying that no energy at all is transferred. 

Early measurements of the accommodation coefficients for rare gases in contact with a 
tungsten surface were carried out by Roberts, although it soon became apparent that his 
experiments were not carried out with sufficiently clean surfaces and thus his data did not 
represent the values expected for the gas-surface interaction with a clean metal.- The work of 
Roberts did, however, stimulate early theoretical investigations, especially for describing the 
interaction of He atoms with surfaces using quantum mechanics.-^ In the 1960s with the 
advent of high vacuum technology and good surface cleaning techniques reliable measure- 
ments of the accommodation coefficient for rare gases on metal surfaces became available 
from two different groups, that of Thomas et al.— ^ and of Kouptsidis and Menzel.— ^ An 
extensive review of work pertaining to accommodation coefficients and a very useful com- 
pendium of available experimental data has been presented by Saxena and Joshi.— Other 
extensive reviews have been given by Goodman and Wachman.— 

The purpose of this paper is to present calculations for the accommodation coefficients 
of the heavy rare gases with metal surfaces using a recently developed classical theory for 
atom-surface collisions that includes both direct scattering and trapping-desorption pro- 
cesses. Similar classical scattering theories have been applied previously to calculations of 
the accommodation coefficient and reasonable agreement with measurements for the heavy 
rare gases on clean tungsten surfaces was obtained, but these calculations included only di- 
rect scattering processes and did not properly include trapping and subsequent desorption by 
the physisorption well of the interaction potential.— The trapping-desorption fraction is that 
portion of an incident beam of gas particles directed towards a surface that gets trapped by 
the physisorption well during the initial collision process. If this fraction remains in the ph- 
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ysisorption well and does not go on to become permanently adsorbed or chemisorbed (which 
is the expected case for rare gas atoms if the temperature is not too small) these physisorbed 
atoms will eventually desorb and the standard assumption is that the trapping-desorption 
fraction leaves the surface in a thermal energy distribution that is nearly in equilibrium at 
the temperature of the surface. Under such an assumption the trapping-desorption fraction 
is expected to enhance the accommodation coefficient and cause it to have values closer to 
unity. On the other hand, the direct scattering fraction tends to exchange less energy with 
the surface and its contribution to the accommodation coefficient is expected to cause it to 
have values less than unity. 

The gas-surface scattering theory applied in this paper uses classical mechanics. In the 
initial collision with the surface a gas particle will either be scattered back into the con- 
tinuum states (the direct scattering fraction) or will be trapped in the physisorption well 
of the interaction potential. Those particles that are trapped can be subdivided into two 
classes, those that have negative total energy and those that have positive total energy 
but are traveling at angles so close to the surface that they cannot escape from the well. 
This latter class is sometimes called the chattering fraction. In the theory used here the 
trapped particles are tracked as they make further collisions with the surface and with each 
subsequent collision a fraction remains trapped but a fraction receives enough energy to 
escape into the continuum states. These subsequent collisions are treated with an iteration 
algorithm that can be carried out to very large numbers until virtually all of the initially 
trapped particles have desorbed. The theory has demonstrated clearly the conditions for the 
trapping-desorption fraction to leave the surface in equilibrium, and it has also explained 
experimental data for rare gas scattering under well-defined conditions for which both a 
direct scattering contribution and a trapping-desorption fraction were observed as separate 
peaks in the energy-resolved scattering spectra.— 

As expected, for the accommodation of He and Ne at a tungsten surface, where quantum 
mechanics should be dominant in the scattering process, the present classical theory is unable 
to explain the measured experimental data. However, good agreement with data is obtained 
for the heavy rare gases Ar, Kr and Xe. 
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II. THEORY 



The energy accommodation coefficient oce{Ts,Tq) is the ratio of the average energy ex- 
changed by a gas in contact with a surface normalized to the maximum thermodynamically 
allowed energy that could be exchanged. 

(T t\ E~f-<E l > E~ f -2k B T G 

MT *' Tg) " <E f >-< Ei > ~ 2k B T s - 2k B T G ■ (1) 

In Eq. ([T]) T G is the temperature of the gas, T s is the temperature of the surface, and k B 
is Boltzmann's constant, Ef is the average energy of a gas particle after making a collision 
with the surface. The expression on the far right hand side of Eq. ([1]) is obtained under the 
assumption that both the gas and surface are in equilibrium at their respective temperatures, 
thus the average energy of the incident gas is < Ei >= 2k B T G and the average energy of 
the gas if it should come into equilibrium with the surface would be < Ef >= 2k B Ts- 
These average energies are obtained from the Knudsen distribution for a gas in equilibrium, 
sometimes called the flux-corrected Maxwell-Boltzmann distribution 
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If the gas is initially in equilibrium then the average final energy after a collision with 
the surface is given by 
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where dR(pf,pi,Ts)/dEf dflf is the differential reflection coefficient giving the probability 
per unit final energy and final solid angle that a gas particle initially in momentum state 
Pi will make a transition to the state p/ as a result of the interaction with the surface. 
The differential reflection coefficient must obey the two conditions of unitarity and detailed 
balancing, as does also the Knudsen distribution of Eq. (J2J). The condition of unitarity 
means that the number of gas particles is conserved, i.e., for a given initial momentum state 
Pi the integral of the differential reflection coefficient over all final energies and angles is 
normalized to unity. 

It is convenient to define an accommodation coefficient that is a function of a single 
temperature by taking the limit as the surface and gas temperatures approach the same 
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value. Thus results in the equilibrium energy accommodation coefficient defined as 

Urn 

oc E {T) = a E (T s ,T G ) . (4) 

T G ^T S ^T 

All calculations in this paper will be for oe(T) since most experimental data for the accom- 
modation of rare gases on clean surfaces is reported in this form. 

Using the condition of detailed balancing, the temperature limit of Eq. (J3J) can be readily 
carried out leading to the final form 



1 r ^ f ,n r^f jn.fv. ^ dP«( Pl ;T G )dR(p f , Pl ;T ; 
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At this point the only remaining quantity needed for evaluating the accommodation 
coefficient is the differential reflection coefficient dR(jpf,pi,Ts)/dEf dQf. This provides 
a complete description of the scattering process which means that it contains not only the 
direct scattering arising from a single collision or a small number of collisions with the surface, 
but it also should contain the contributions of those particles that are initially trapped and 
then subsequently desorbed. The present authors have recently developed a complete theory 
of atom-surface scattering that includes both contributions.— This theory is based on an 
initial differential reflection coefficient for a single surface collision dR (p/, Pi, T$)/dEf dQf. 
This initial collision results in a scattered intensity that consists of a fraction that is the 
direct scattering contribution which leaves the surface and a fraction that is trapped in the 
physisorption well of the interaction potential. The trapped fraction is followed inside the 
well and when those particles have another collision, some escape into the continuum and 
some remain trapped and go on to have further collisions with similar consequences. By 
dividing all trapped particles into a distribution of small energy and angular bins they can 
be followed through many collisions until ultimately essentially all of them have escaped into 
the continuum and desorbed. This process can be written schematically as the following 
equation 

dR(p f ,pi) dR°(p f ,pi) r dR°(p fl p b ) dR°(p b ,pj) 

dE f dft f ~ dEfdflf J b b dE f dtt f dE b dQ b 1 ' 

+ fdE b dn b dR °^ Pb) dRl b»>*) + ... 
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where dR n (pb, p i: Ts)/dEb dflb is the differential reflection coefficient giving the distribution 
of particles remaining trapped in the bound states after n collisions and the intermediate 
integrations in the higher order terms are carried out only over angles and energies that 
pertain to the trapped fraction. The process described in Eq. (jH]) is readily developed into 
an iterative scheme in which the scattered distribution remaining in the well after the nth 
collision becomes the source for the next collision. The details of the calculation of the 
differential reflection coefficient of Eq. are given in Ref.p^] where it is shown that this 
procedure for calculating the trapping-desorption fraction can not only explain the physical 
behavior of the trapping-desorption fraction, but it can also explain experimental data for 
energy-resolved scattering spectra of rare gases taken under conditions in which there are 
distinct and clearly exhibited features due to both direct scattering and trapping-desorption. 

The zeroth order differential reflection coefficient is chosen to be that due to a potential 
consisting of an attractive square well in front of a smooth repulsive wall whose surface vi- 
brates due to the thermal motion of the substrate atoms. The square well has depth D and 
width b. For a classical mechanical treatment a square well is a reasonable approximation 
since it describes correctly the increase in energy and refraction towards more normal inci- 
dence angles when a particle enters the physisorption potential. The width b does not affect 
the scattered intensities provided it is larger than the selvedge region of the surface, i.e., 
as long as it is larger than the surface vibrational displacements. However, trapping times 
are proportional to b. This differential reflection coefficient has been shown to be useful in 
explaining a variety of atom- surf ace scattering experiments and is given by2&£L22£LM. : 

dR°{p f , Pi ;T s ) 
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where AEq = (p/ — Pi) 2 /2M is the recoil energy. pi Z is the z component of the incident 
momentum, |t/j| 2 is the form factor of the scattering is the normalization coefficient, 
P = (P/ — Pj) is the parallel momentum exchange, and S u . c , is the area of a surface unit cell. 
The quantity Vr is a weighted velocity of sound parallel to the surface. It is expected to have 
a values of order of the Rayleigh phonon speed and can be calculated from the complete 
surface phonon spectral densiy, however, it is usually taken to be a constant.— ■2^^- 

The amplitude r/j of the form factor appearing in Eq. ([7]) is in general the transition 
matrix of the elastic interaction potential extended off of the energy shell to account for 



inelastic scattering. A good approximation that has been extensively used is the first-order 
distorted wave Born approximation matrix element which for a strongly repulsive surface is 
given by 17 

r fi = Ap fz p iz /m (8) 

The main numerical operations involved in carrying out calculations are the multiple 
integrals involved in the accommodation coefficient of Eq. (j3J) and in the iterative evaluation 
of the differential reflection coefficient of Eq. ([6]). In each case these are six-dimensional 
integrations, although, if the surface is azimuthally symmetric as is the case for the potential 
used here, the accommodation coefficient reduces to a five- dimensional integral. The angular 
integrations are carried out using Gauss- Legendre algorithms and the energy integrals with 
Gauss-Laguerre algorithms. 

III. RESULTS 

Comparisons with experimental data for calculations using the theory and interaction 
potential described above are presented for the heavy rare gases in contact with a tungsten 
surface in Figures [Ul3l The data from Thomas et al. are shown as open circles and the data 
from Kouptsidis and Menzel are shown as filled circles. 

Fig. [1] shows the measured equilibrium accommodation coefficient for Ar on W compared 
to four curves calculated with different well depths of 10, 20, 25 and 50 meV. The velocity 
parameter is Vr = 500 m/s. The best agreement with the data is for a well depth of ap- 
proximately D = 25 meV. A table of measured and theoretically calculated well depths for 
the Ar/W system is given in Ref. p 8 -] which shows that this value of 25 meV is somewhat 
smaller than expected. This table is based on values presented in Ref. p&] and the measure- 
ments, primarily obtained from thermal desorption experiments, range from 78 to 127 meV 
while calculated values are somewhat smaller ranging from 33 to 47 meV. If the velocity 
parameter vr is made somewhat larger the calculations for a larger well depth approach 
more closely to the data at large temperatures, but the agreement at low temperatures be- 
comes worse. Although the well depth predicted by these calculations is somewhat smaller 
than expected, it is considerably larger than the value of 15 meV used previously to fit the 
data with calculations based solely on the direct scattering. 18 Thus it becomes clear that 
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including the trapping-desorption in the calculation significantly increases the value of the 
accommodation. 

Fig. [2] shows the accommodation coefficient data for Kr/W compared to calculations 
carried out for two different well depths, 50 and 70 meV. The best agreement with the data 
is for a well depth of approximately D = 50 meV. Larger well depths lead to larger trapping- 
desorption fractions, and since the trapping-desorption fraction is nearly in equilibrium this 
tends to enhance the accommodation coefficient. As in the case for Ar/W, calculations 
with vr larger than 500 m/s will tend to decrease the accommodation coefficient for a given 
well depth, but at the expense of poorer overall agreement with the data. Estimated well 
depths for the Kr/W system have been obtained only from thermal desorption experiments 
and these values range from 195 to 247 meV as tabulated in Ref. Thus the value 

used here to give a best fit with the data is small in comparison to the thermal desorption 
measurements, but again as for the Ar/W system it is twice as large as the value obtained 
for calculations based only on direct scattering.— 

Fig. [3] shows similar comparisons with data for the case of Xe/W. Calculations for two 
well depths, 100 and 150 meV, are shown and Vr = 500 m/s. Both of these well depths are 
somewhat smaller than the independent measured thermal desorption value of 180 meV. 
It is interesting to note that the calculations for D = 150 meV and temperatures below 
Ts = 150 K show that essentially all of the gas atoms are trapped in the physisorption well 
and escape nearly in equilibrium which results in complete accommodation, or = 1. 

It is to be expected that a classical theory will not be adequate to describe the accom- 
modation of the lighter rare gases He and Ne on a surface of atoms as heavy as tungsten. 
Numerous treatments have shown that the interaction of these gases with metal surfaces, 
especially for the case of He gas, is quantum mechanical in nature and the scattering is dom- 
inated by elastic and single phonon inelastic processes. A classical theory, such as used here, 
cannot properly treat quantum mechanical processes and the present calculations predict 
accommodation coefficient values that are much too large for He and Ne. However, there are 
several quantum theoretical treatments of the accommodation coefficient based on exchange 
of small numbers of phonons which explain the measured values for the He/W and Ne/W 
system quite nicely.— 1 ^ 
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IV. CONCLUSIONS 



This paper presents calculations of the equilibrium accommodation coefficient for energy 
exchange at a gas-surface interface using a newly developed theory of atom-surface scattering 
in the classical limit that treats both the direct scattering and the fraction of particles that 
are trapped and subsequently desorbed after the initial collision with the surface. This theory 
is applied to a relatively straightforward model of the interaction potential, consisting of a 
strongly repulsive vibrating repulsive wall with an attractive square physisorption well in 
front. However, the theory treats the statistical mechanics of the scattering process properly 
and is able to track all initially physisorbed particles until they eventually desorb. This 
theory not only describes correctly both the direct and trapping-desorption fractions, it has 
been used to explain measured experimental data whose energy-resolved scattering spectra 
exhibit distinct features due to direct and trapping-desorption events. Thus, it is of interest 
to calculate the accommodation coefficient using this theory to see if it explains the available 
data for energy transfer at a gas-surface interface. 

A large amount of data exists for the accommodation of a variety of atomic and molecular 
gases at different types of surfaces. However, the most carefully defined systems, both 
experimentally and theoretically, are the rare gases accommodating at a tungsten surface. 
Although data is available for all the rare gases except radon, comparisons here are made 
only for the heavier rare gases. This is because the light mass rare gases, He and Ne, interact 
quantum mechanically and are not well explained by a purely classical theory. 

This work can be viewed as a logical extension of an earlier paper by one of the authors 
in which calculations with a similar interaction potential model, but a theory that contained 
only the direct scattering component, was applied to the energy accommodation coefficient.— 
Thus, the present work when compared to the previous results gives a clear indication of 
the contributions of the trapping-desorption fraction to the accommodation. 

Good overall agreement between calculations and measured accommodation coefficient 
data is obtained. However, the results do depend on the choice of the well depth and the 
velocity parameter that arises from the model of the interaction potential. Neither of these 
quantities has been well established for the interaction of heavy rare gases with the tungsten 
surface. The calculated values of the well depths that give the best agreement with mea- 
surements tend to be somewhat smaller than estimates extracted from thermal desorption 
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experiments, although there are typically significant difference between such measurements 
in the cases where more than one exists. In comparison with the previous calculations, how- 
ever, the present work predicts well depths that are significantly larger due to the influence 
of the trapping-desorption fraction. 

The fact that this theory explains the available data for heavy rare gas accommodating 
at clean tungsten surfaces, and the fact that state-to-state calculations explain recently 
available data for Ar scattering under conditions where the energy-resolved spectra exhibited 
clear evidence for distinct direct scattering and trapping-desorption features implies that it 
should be useful for predicting the energy accommodation for other gas-surface systems. In 
particular it should be able to predict the behavior of other systems as a function of the 
experimentally accessible initial conditions such as temperature, well depth, gas particle 
mass and surface mass. 
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FIG. 1: The equilibrium energy accommodation coefficient o.e{T) as a function of the temperature 
Ts for Ar on a W surface. Calculations four with different well depths as marked are compared 
with the experimental data of Thomas et a/ i 10 ' 11 (open circles) and Kouptsidis and Menze l 12 ' 13 
(filled circles). 



13 



o 
< 

LU 



1 



0.8 



0.6 



0.4 



0.2 



EAC vs. Temperature (Continuum) 

Kr/W; v R =500mps 

1 1 1 1 1 

N 

\ 

\ 

\ 

O \ N 



'5, 



O 



D=50meV 
D=70meV 
o Data 



o 



o 



100 200 
T S (K) 



300 



FIG. 2: cce(T) as a function of the temperature Tg for krypton gas in contact with a tungsten 
surface. Two calculations with well depths of 50 and 70 meV are shown as marked and the data 
points are from the same sources and are denoted the same as in Fig. [TJ 
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FIG. 3: The equilibrium energy accommodation coefficient o.e{T) as a function of the temperature 
Tg for xenon in contact with a tungsten surface. Data are denoted as in Fig. [TJ Calculations with 
well depths 100 and 150 meV are shown. 
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